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Abstract 

Detailed  comparisons  of  fluid  flow,  concentration,  and  electric  fields  between  the  planar  solid  oxide  fuel  cell  (planar  SOFC)  and  the  mono¬ 
block-layer-built  solid  oxide  fuel  cell  (MOLB-type  SOFC)  are  made  in  the  present  study.  A  multi-physics  model  is  developed  to  simulate 
the  SOFCs  of  different  geometries.  The  model  couples  the  electrochemical  kinetics  with  the  fluid  dynamics  and  multi-component  species 
transport.  Darcy’s  law  governs  the  flow  momentum  in  the  porous  electrodes.  The  Bruggman  model  depicts  the  diffusion  of  reactant  species 
within  the  electrodes.  The  electrochemical  kinetics  and  the  species  transport  on  the  electrode  surfaces  is  connected  by  the  Butler- Volmer 
equation.  Results  show  that  the  MOLB-type  SOFC  has  a  higher  fuel/oxidant  utilization  than  the  planar  SOFC.  The  high  utilization  of  the 
fuel/oxidant  reflects  the  high  current  through  the  SOFC. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC)  is  a  potential  candi¬ 
date  for  future  energy  conversion  schemes  because  of  its 
higher  conversion  efficiency  than  the  conventional  heat  en¬ 
gine.  Besides,  internal  reforming  of  hydrocarbon  fuels  can 
be  performed  at  the  anode,  which  allows  direct  use  of  nat¬ 
ural  gas,  LPG,  and  methanol  as  fuels.  Underlying  such  pro¬ 
jected  benefits  are  its  fuel  flexibility  and  higher  efficiency 
compared  to  the  conventional  energy  conversion  system  such 
as  ICEs. 

There  is  a  large  amount  of  research  on  SOFCs.  Fo¬ 
cus  has  been  placed  on  the  designs  of  new  material  [1-4] 
and/or  novel  geometry  [5-8].  The  planar  SOFC  has  received 
much  more  attention  than  the  tubular  SOFC  because  of 
its  easier  fabrication  and  higher  power  density.  Fig.  1(a) 
schematically  shows  a  typical  planar  SOFC.  PEN  (posi¬ 
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tive  electrode-electrolyte-negative  electrode)  plate  is  sand¬ 
wiched  between  two  inter-connectors  (bipolar  plates).  The 
inter-connectors  are  machined  with  passages  for  channeling 
the  fuel  and  air  to  the  fuel  cell.  Fig.  1(b)  shows  another  ge¬ 
ometry  known  as  the  MOLB  (mono-block-layer-built)-type 
SOFC.  The  PEN  plate  is  molded  into  a  corrugated  shape.  The 
fuel  and  air  run  in  the  trapezoidal  channels  that  are  separated 
by  the  corrugated  PEN  plate.  Since  the  inter-connectors  are 
flat,  the  laborious  work  of  channel  machining  is  not  neces¬ 
sary. 

The  objective  of  the  present  study  is  to  compare  the  trans¬ 
port  and  electrochemical  characteristics  between  the  conven¬ 
tional  planar  SOFC  and  the  MOLB-type  SOFC.  A  CFD- 
based,  multi-physics,  and  electrochemical-transport  coupled 
SOFC  model  is  developed.  The  model  involves  solving  con¬ 
servation  equations  for  mass,  momentum,  energy,  multi- 
component  species,  and  electrical  charge  along  with  elec¬ 
trochemical  kinetics  both  in  the  active  and  inactive  regions 
of  the  SOFC.  The  results  obtained  by  this  study  can  assist  the 
designer  in  understanding  how  geometry  affects  gas  trans¬ 
port  and  electrochemical  characteristics  in  the  SOFC.  Also, 
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Nomenclature 

a ' 

stoichiometric  coefficient  of  the  products 

a " 

stoichiometric  coefficient  of  the  reactants 

D 

diffusivity  [m2  s-1] 

F 

Faraday’s  constant  [96487  C  mol-1] 

• 

1 

current  density  [Am-2] 

• 

J 

transfer  current  density  [A  m-3] 

j 

diffusive  flux  [kg  s-1  m-2] 

k 

thermal  conductivity  [WK-1  m-1] 

M 

molecular  weight  [kg  mol-1] 

P 

pressure  [Pa] 

Pi 

possibility  of  the  same  particles  stretch  from 

the  inter-connector  to  the  electrolyte 

R 

universal  gas  constant  [W  mol-1  K-1] 

T 

temperature  [K] 

U 

velocity  vectors  [ms-1] 

X,  Y,  Z 

coordinate  system,  Fig.  1  [m] 

Xi 

fraction  of  molar  concentration  of  species  i 

Yi 

mass  fraction  of  species  of  i 

Greek  symbols 

a 

symmetric  factor 

£ 

porosity 

r)  (or  Over_P)  overpotential,  @c  ~  [V] 

K 

permeability  [m2] 

[i 

viscosity  [ms-2] 

P 

density  [kgm3] 

0 

electric  conductivity  [£2-1  m-1] 

r 

tortuosity  of  the  porous  electrode 

0 

solid-phase  potential  [V] 

Subscripts 

A 

anode 

C 

catalyst  or  cathode 

eff 

effective 

E 

electrolyte 

F 

fluid 

i 

species 

ref 

reference 

S 

solid 

T 

transfer  current 

they  can  provide  a  more  complete  basis  for  optimizing  the 
geometry  of  the  SOFC  stack. 

2.  Modelling 

2.1.  Basics  of  SOFC 

The  oxygen  reduction  reaction  on  the  porous  cathode 
(LSM/YSZ)  of  the  SOFC  is 


The  oxygen  ion  moves  to  the  anode  (nickel/YSZ)  through  the 
electrolyte  (YSZ).  It  then  reacts  with  either  H2  or  CO  on  Ni 
catalysts,  i.e. : 

72— 


h2  +  o 
co  +  o; 


H2O  +  2e“ 
C02  +  2e' 


(2) 

(3) 


The  co-oxidations  of  Eqs.  (2)  and  (3)  are  currently  limited  to 
a  single-step  reaction,  i.e.: 

72— 


H2  +  CO  +  20' 


H20  +  C02  +  4e 


(4) 


2.2.  Model  equations 


The  PEN  plate  is  a  porous  medium,  while  the  flow  chan¬ 
nels  are  filled  with  reactant  gases  only.  In  this  model,  the 
entire  domain  is  treated  as  a  porous  matrix  with  proper  porosi¬ 
ties  assigned  to  each  sub-domain.  For  example,  the  porosity 
of  the  open  channel  is  set  to  1,  and  extremely  small  values  of 
porosity  are  assigned  to  the  electrolyte  and  inter-connectors. 
As  a  result,  the  entire  domain  can  be  governed  by  a  single  set 
of  equations.  Thus,  it  is  not  necessary  to  specify  the  boundary 
conditions  at  the  interfaces  between  sub-domains. 

2.2.1.  Gas-phase  transport  and  reaction 

The  mass  and  momentum  conservations  are  described  as 
the  following  equations: 


V  •  (epV)  =  0 


V  •  (epUU)  =  —  eV p  +  V  •  (£/xeffVU)  + 


£2p\] 


K 


(5) 


(6) 


As  for  the  energy  conservation,  it  is  assumed  that  the  fluid 
and  solid  phases  are  in  local  thermal-equilibrium,  and  thus 
the  thermal  dispersion  is  negligible.  In  addition,  the  thermo¬ 
physical  properties  of  both  solid  and  fluid  are  constant.  The 
equation  governing  the  energy  conservation  is  expressed  as 


V  •  (spCpllT)  —  V  •  (keffVT)  —  j jij  + 


1  •  1 


^cff 


(7) 


In  the  above  equation,  the  radiation  heat  flux  and  energy  dissi¬ 
pation  due  to  the  shear  stress  are  neglected.  The  effective  ther¬ 
mal  conductivity  can  be  determined  by  the  following  equation 

[9]: 


keff  —  —  2ks  + 


1 


(e/(2*s +**))  + ((1  -  e)/3*s) 


(8) 


where  kp  is  the  fluid-phase  conductivity,  and  ks  the  weight- 
averaged  conductivity  between  the  ionic  conductor  and  the 
electric  conductor.  The  last  two  terms  of  Eq.  (7)  represent  the 
effects  of  conversion  loss  and  joule  heating,  respectively. 

The  conservation  of  species  is  expressed  as 


V  •  (epUYi)  =  V  •  J;  +  (a'{  -  a^ 

r 


/  Jt 


(9) 


O2  +  4e 


20' 


(1) 


a"  and  o'  are  the  stoichiometric  coefficients  of  the  products 
and  reactants,  respectively.  The  diffusive  mass  flux  of  the 
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(a)  planar  SOFC 


PEN 

plate 


Bipolar  plate 

Cathode  (LSM/YSZ) 
Electrolyte  (YSZ) 
Anode  (Ni/YSZ) 

Bipolar  plate 


(b)  MOLB-type  SOFC 


PEN 

plate 


Bipolar  plate 


Cathode  (LSM/YSZ) 
Electrolyte  (YSZ) 

Anode  (Ni/YSZ) 

Bipolar  plate 


Fig.  1.  Schematic  drawing  of  the  planar  and  MOLB-type  SOFCs. 


species  J;  is  modeled  as 

J,  =  pA-.effVF;  +  ^  A.effVM  -  M  J^Dj^Yj 

j 

-vmJ2  djwYj  (10) 

i 

The  first  term  on  the  right-hand  side  represents  the  Fickian 
diffusion  due  to  concentration  gradient.  The  last  three  terms 


are  the  corrections  to  enforce  the  Stefan-Maxwell  equations 
for  the  multi-component  diffusion  system.  The  effective  dif- 
fusivities  in  a  multi-component  mixture  of  the  porous  layers 
follow  the  Bruggman  model  [10]: 

DiMf  =  8TDi  (11) 

The  relationship  between  the  electrochemistry  and  the 
species  transport  kinetics  on  the  electrode  surfaces  is  depicted 
by  the  Butler- Volmer  equation  [11].  That  is,  the  transfer  cur¬ 


ia)  Anode  (b)  Cathode 


Fig.  2.  Current  transports  in  the  electrode  of  the  SOFC. 
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rent  densities  in  the  anode  and  cathode,  respectively,  are: 


Jt,a  —  Jo, a  ' 


1/4 


^H2,ref 


/  Xn2p 

\^H20,ref 


1/4 


/  *CO  \ 
\^CO,ref/ 

/  *C02 
\^C02,ref, 


1/4 


exp 


■1/4 


oaF 

RT 


0 


x  exp 


-(l-aA)F 
RT 


0 


(12) 


Jt,c  —  Jo,c 


Xo, 


1/4 


^02,ref . 


exp 


-etc  F 

RT 


T] 


(13) 


where  jo  is  the  exchange  current  density,  r\  the  electrode  over¬ 
potential,  a  a  and  etc  the  anode  and  cathode  Tafel  constants, 
respectively. 


2.2.2.  Solid-phase  conduction 

The  porous  electrode  of  the  SOFC  contains  two  kinds  of 
solid-phase  material,  i.e.,  electrolyte  (ionic  conductor,  e.g. 


Table  1 


Geometric  comparison  between  the  planar  and  MOLB-type  SOFCs 


Component 

Parameters 

Planar 

SOFC 

MOLB-type 

SOFC 

Anode 

Thickness  (mm) 

1.4 

1.4 

Length  (mm) 

12 

16.8 

Electrolyte 

Thickness  (mm) 

0.2 

0.2 

Length  (mm) 

12 

16.8 

Cathode 

Thickness  (mm) 

0.4 

0.4 

Length  (mm) 

12 

16.8 

Fuel  channel 

Cross-sectional  area  (mm2) 

90.24 

89.6 

Length  (mm) 

100 

100 

Air  channel 

Cross-sectional  area  (mm2) 

90.24 

89.6 

Length  (mm) 

100 

100 

Module 

Volume  (mm3) 

28000 

28300 

YSZ)  and  catalyst  (electric  conductor,  e.g.,  LSM  or  Ni).  A 
potential  difference  exists  between  the  catalyst  (<Pq)  and  elec¬ 
trolyte  (0e)  to  drive  the  transfer  current  (jx),  keeping  the 
electrochemical  reaction  continuously.  As  shown  in  Fig.  2, 


Table  2 


Porous-electrochemical  properties  of  the  electrodes 


Description 

Unit 

Value 

Exchange  current 

Anode,  j0,A 

Am-3 

1.0  X  109 

Cathode,  j0,c 

Am-3 

5.0  x  107 

Transfer  coefficient 

Anode,  a  a 

— 

1.0 

Cathode,  ac 

— 

1.0 

Porosity,  e 

Anode 

— 

0.4 

Cathode 

— 

0.5 

Electrolyte 

— 

1.0  X  10“5 

Permeability,  k 

Anode 

7 

m 

1.0  x  10“12 

Cathode 

2 

m 

1.0  x  10“12 

Electrolyte 

2 

m 

1.0  x  10“18 

Tortuosity,  t 

Anode 

— 

1.5 

Cathode 

— 

1.5 

Electrolyte 

— 

6.0 

Effective  electric  conductivity,  oreff 

Anode 

Q~[  m-1 

1.0  X  105 

Cathode 

f}-1  m-1 

7.7  x  103 

Electrolyte 

Q~l  m-1 

1.0  x  10“2° 

Inlet  pressure 

Anode 

kPa 

1.013  x  105 

Cathode 

kPa 

1.013  x  105 

Stack  temperature 

Anode 

K 

973 

Cathode 

K 

973 

Mass  flow  rate 

Anode 

kgs-1 

1.19  x  10“6 

Cathode 

kgs-1 

2.98  x  10“5 

Species  molar  concentration  at  anode  inlet 

Hydrogen,  Xr2 

— 

45.57% 

Carbon  monoxide,  Xqo 

— 

6.97% 

Carbon  dioxide,  Xqo2 

— 

6.55% 

Water,  Vh9o 

— 

36.84% 

Methane,  Xch4 

— 

4.07% 

Total 

— 

100% 

Species  molar  concentration  at  cathode  inlet 

Oxygen,  X0 

— 

21.11% 

Nitrogen,  X n 

— 

78.99% 

Total 

— 

100% 
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Fig.  3.  Comparison  of  the  secondary-flow  patterns  between  the  planar  and  MOLB-type  SOFC. 


the  current  passes  through  electrode  can  be  decomposed  two 
parts,  i.e. : 

i  =  ic  +  iE  (14) 

ic  and  ig  are  the  currents  flowing  through  the  catalyst  and  the 
electrolyte,  respectively.  Since  the  electrodes  are  electroneu¬ 
tral  everywhere,  there  is  no  charge-buildup  in  the  active  lay¬ 
ers.  Thus,  the  charge  conservation  is 

V  •  i  =  0  (15) 

That  is 

V  •  ic  =  —V  •  iE  (16) 

These  two  current  components  interact  through  electrochem¬ 
ical  reactions.  The  electrons  are  transferred  to  the  catalyst 
from  the  electrolyte  in  the  anode,  and  vice  versa  in  the  cath¬ 


ode.  Application  of  Ohm’s  law  to  Eq.  (16)  yields  the  current 
conservation: 

V  •  (—cre ff,cV^c)  =  -J*t  (I7) 

v  •  (-^efFE^^E)  =  jT  (18) 

where  <reff,c  and  <xcfFE  are  the  effective  electric  conductivities 
of  the  catalyst  and  electrolyte,  respectively.  They  are  modeled 
as 


=  erC(l  -  s)<pcPi,C 

(19) 

Ceff.E  =  OE(l  —  £)<pEPi,E 

(20) 

where  0c  and  0e  are  the  volume  fraction  of  the  catalyst  and 
electrolyte,  respectively. pt£  is  the  possibility  of  the  catalyst 
in  the  connection  of  the  inter-connector,  and  is  the  pos¬ 
sibility  of  electrolyte  in  the  connection  of  PEN  plate  [12].  It 
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Fig.  4.  Mass  fraction  distributions  of  the  oxygen  and  hydrogen  in  the  SOFCs:  (a)  planar  SOFC;  (b)  MOLB-type  SOFC. 


is  noted  that  only  a  long-range  connection  of  the  same  parti¬ 
cles  stretch  from  the  inter-connector  to  the  electrolyte  ensures 
good  conductivity. 

2.2.3.  Geometric  parameters  and  physical  properties 

Table  1  lists  the  geometric  parameters  of  the  planar  SOFC 
and  the  MOLB-type  SOFC  in  a  single  module.  To  place  a 
common  basis,  several  important  parameters  are  kept  at  the 
same  value  between  these  two  geometries  such  as  the  PEN- 
plate  thickness,  the  channel  cross-sectional  area  and  length, 
and  the  module  volume.  The  majority  of  discrepancy  between 
these  two  designs  is  the  length  of  the  PEN  plate.  The  MOLB- 
type  SOFC  has  a  16.8  mm  PEN  plate  in  length,  which  is  about 
40%  longer  than  that  of  the  planar  SOFC. 

Other  parameters  and  physical  properties  used  in  this 
model  are  summarized  in  Table  2.  The  ambient  air  supplies 
the  cathode  with  oxidants.  The  anode  is  fed  with  a  simu¬ 
lated  reformed  gas  i.e.  a  composition  of  CH4,  CO,  CO2, 
N2,  and  H2O.  The  molar  concentration  of  the  anode  fuel 


is  45.57/6.97/6.55/36.84/4.07%  for  H2/CO/CO2/H2O/CH4. 
Among  them,  CH4  and  N2  are  considered  as  inert  gases  and 
serve  as  diluents.  The  operation  voltage  of  the  cell  is  set  to 
0.55  V,  meaning  that  the  overpotential  within  the  module  is 
0.4  V  with  an  OCV  of  0.95  V.  In  addition,  the  fluid  viscosity 
(/reff)  is  calculated  using  mixture  kinetic  theory,  specific  heat 
(Cp)  via  JANNAF  tables  for  mixtures,  and  thermal  conductiv¬ 
ity  (kp)  is  calculated  for  constant  Prandtl  number  (Pr  =  0.7). 

3.  Numerical  implementation 

The  governing  equations  are  numerically  solved  by  the 
control- volume-based  finite  difference  method  [13].  The 
momentum  equations  are  solved,  followed  by  a  pressure 
correction  equation  that  does  the  mass  balance.  Species 
transport  equations  and  electric  conduction  equations  are 
solved  simultaneously  after  the  bulk  flow  calculation.  Veloc¬ 
ity  control  volumes  are  staggered  with  respect  to  the  main 
control  volumes,  and  coupling  of  the  pressure  and  velocity 
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Fig.  5.  Current  density  distributions  along  the  solid-phase  of  SOFCs:  (a)  planar  SOFC;  (b)  MOLB-type  SOFC. 


fields  is  treated  via  the  SIMPLER  pressure  correction  algo¬ 
rithm  [14-16].  Solutions  are  considered  to  be  converged  at 
each  test  condition  after  the  ratio  of  residual  source  (includ¬ 
ing  mass,  momentum,  species,  and  charge)  to  the  maximum 
flux  across  a  control  surface  becomes  below  1.0  x  10~6. 
Computations  are  performed  on  structured  grids  with  number 
of  12  x  44  x  50  and  18  x  40  x  50  for  the  planar  SOFC  and 
the  MOLB-type  SOFC,  respectively.  Grid-independence 
test  is  carried  out  by  running  three  additional  coarser  and 
finer  meshes.  A  typical  simulation  requires  about  300  min  of 
central  processing  unit  time  on  a  Pentium  IV  2.8  GHz  PC. 

4.  Results  and  discussion 

4.1.  Secondary  flow 

Fig.  3  compares  the  secondary-flow  pattern  in  the  flow 
channel  and  the  porous  electrode  between  the  planar  SOFC 


and  the  MOLB-type  SOFC.  The  cross-sectional  view  is  cut 
through  the  mid-plane  of  the  computational  module.  Al¬ 
though  the  structure  of  the  secondary  flow  between  these 
two  geometries  is  different,  the  general  trend  in  the  move¬ 
ment  of  the  secondary  flow  is  similar  in  essential.  It  is  briefly 
described  as  follows. 

In  the  anode,  the  secondary  flow  is  directed  from  the  elec¬ 
trolyte  toward  the  flow  channel.  This  is  because  the  gas  bulk 
density  in  the  anode  is  increased  due  to  the  electrochemi¬ 
cal  oxidation  reaction  (Eq.  (4)).  Thus,  the  outflow  should  be 
higher  than  the  inflow  in  the  anode  for  mass  conservation. 
In  contrast,  the  oxygen  reduction  reaction  in  the  cathode  re¬ 
duces  the  air  bulk  density  (Eq.  (1)).  Thus,  the  fresh  air  flows 
continuously  into  the  cathode  from  the  channel  to  satisfy  the 
continuity.  In  the  channel  flow,  the  velocity  vectors  of  the 
secondary  flow  are  always  directed  toward  the  channel  cen¬ 
ter  due  to  the  wall  effect.  The  stronger  secondary  flow  in  the 
cathode  channel  is  because  of  the  higher  inlet  velocity  of  the 
oxidant. 
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4.2.  Species  concentration 

A  comparison  of  the  distributions  of  Yo2  (oxygen)  and 
Yr2  (hydrogen)  between  the  planar  SOFC  and  the  MOLB- 
type  SOFC  is  shown  in  Fig.  4.  The  general  trend  of  the  con¬ 
centration  distributions  for  both  types  of  SOFCs  is  briefly 
described  as  follows.  Since  the  anode/cathode  acts  a  species 
sink  to  descend  the  fuel  gas/oxidant,  the  gradient  of  species 
concentration  in  the  both  electrodes  is  significant.  In  the  flow 
channel,  however,  the  species  concentration  is  relatively  uni¬ 
form.  In  the  cathode  of  the  planar  SOFC,  a  significant  gra¬ 
dient  of  Yo2  is  found  in  the  lateral  (X)  direction.  The  local 
minimum  of  Yo2  is  found  in  the  region  in  contact  with  the 
inter-connector.  This  is  because  the  fresh  air  cannot  access 
to  this  region  easily.  In  the  anode  of  the  planar  SOFC,  the 
gradient  of  7h?  distribution  is  significant  along  the  normal 
(T)  direction  but  no  so  significant  along  the  lateral  direction. 

As  for  the  MOLB-type  SOFC,  the  local  minimum  of  To2 
is  located  at  the  middle  of  the  upper  plateau  in  contact  with 
the  inter-connecter.  The  local  minimum  of  Th2  occurs  at  the 
corner  formed  by  the  inclined  plane  and  the  ceiling  plane 
where  the  hydrogen  is  transported  mainly  through  diffusion. 

It  is  seen  form  Fig.  4  that  the  mean  concentration  of  H2/O2 
in  the  channel  flow  of  the  MOLB-type  SOFC  is  lower  than 
that  of  the  planar  SOFC.  This  means  that  the  fuel/oxidant 
utilization  is  higher  in  the  MOLB-type  SOFC  than  that  in 
the  planar  SOFC.  The  higher  fuel/oxidant  utilization  is  be¬ 
cause  the  longer  PEN  plate  is  accompanied  by  the  MOLB- 
type  SOFC.  The  longer  PEN  plate  provides  a  more  active 
surface  for  the  electrochemical  reaction.  The  significant  fuel 
utilization  will  reflect  the  high  current  through  the  fuel  cell. 

4.3.  Electrochemical  characteristics 
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Fig.  5(a)  and  (b)  shows  the  distributions  of  the  current 
density  along  the  solid  phase  of  the  planar  SOFC  and  the 
MOLB-type  SOFC,  respectively.  The  current  is  driven  from 
the  anode  to  cathode  inside  the  fuel  cell.  The  arrows  indicate 
the  direction  as  well  as  the  magnitude  of  the  currents.  Quali¬ 
tatively,  these  two  plots  have  several  similar  points.  First,  the 
current  density  gradually  decreases  along  its  direction  in  the 
anode  due  to  the  electrochemical  oxidation  reaction  (Eq.  (4)). 
Conversely,  the  current  density  gradually  increases  along  its 
direction  in  the  cathode  due  to  the  current  accumulation  by 
oxygen  reduction  reaction  (Eq.  (1)).  Secondly,  the  anode  has 
a  lower  current  density  than  the  cathode  due  to  its  higher  area 
of  the  current  passage.  The  major  difference  between  these 
two-type  SOFCs  is  the  magnitude  of  the  current  density.  The 
MOLB-type  SOFC  has  higher  current  densities  than  the  pla¬ 
nar  SOFC. 

Fig.  6  compares  the  overpotential  distributions  on  the  ac¬ 
tive  regions  between  the  planar  SOFC  and  the  MOLB-type 
SOFC.  In  the  anode,  the  positive  overpotential  directs  the 
transfer  current  from  the  catalyst  to  the  electrolyte  (Fig.  2(a)), 
and  vice  versa  in  the  cathode  (Fig.  2(b)).  In  general,  the  cath¬ 
ode  requires  a  higher  energy  to  activate  the  reactants  than 


Fig.  6.  Overpotential  distributions  on  the  surfaces  on  the  solid  phase  of 
SOFCs:  (a)  planar  SOFC;  (b)  MOLB-type  SOFC. 

the  anode.  Therefore  the  absolute  value  of  the  overpotential 
in  the  cathode  is  much  higher  than  that  in  the  anode.  In  ad¬ 
dition,  it  is  higher  in  the  region  of  lower  hydrogen/oxygen 
concentration  (Fig.  5).  That  is  depleting  the  reactants  results 
in  increasing  the  overpotential. 

5.  Conclusions 

The  present  study  has  made  a  detailed  comparison  of  the 
transport  and  electrochemical  characteristics  between  the  pla¬ 
nar  SOFC  and  the  MOLB-type  SOFC.  A  complete  set  of  con¬ 
servation  equations  for  mass,  momentum,  species,  energy, 
and  charge  are  numerically  solved  with  proper  account  of 
electrochemical  kinetics.  The  model  implements  the  voltage- 
to-current  algorithm,  coupling  the  potential  field  with  the 
species  concentration  field.  A  realistic  spatial  variation  of 
electrochemical  kinetics  is  obtained  by  simultaneously  solv¬ 
ing  the  charge  transport  equations  of  the  catalyst  and  elec¬ 
trolyte  in  the  electrode  for  the  first  time.  Informative  results 
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of  the  reactant  gas  concentration  distribution,  overpotential 
distribution,  and  current  density  distribution  have  been  pre¬ 
sented  and  compared.  It  is  found  that  the  MOLB-type  SOFC 
has  a  higher  fuel/oxidant  utilization  than  the  planar  SOFC 
due  to  its  longer  PEN  plate.  The  significant  utilization  of  the 
fuel/oxidant  reflects  the  high  current  passing  through  the  fuel 
cell. 
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